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We develop a unified scattering approach to dynamical Casimir problems which can be applied 
to both accelerating boundaries, as well as dispersive objects in relative motion. A general (trace) 
formula is derived for the radiation from accelerating boundaries. Applications are provided for 
objects with different shapes in various dimensions, and undergoing rotational or linear motion. 
Within this framework, photon generation is discussed in the context of a modulated optical mirror. 
For dispersive objects, we find general results solely in terms of the scattering matrix. Specifically, 
we discuss the vacuum friction on a rotating object, and the friction on an atom moving parallel to 
a surface. 



I. INTRODUCTION 

Quantum zero-point fluctuations manifest themselves 
in a variety of macroscopic effects. A prominent example 
is Casimir's demonstration that these fluctuations lead 
to attraction of two perfectly conducting parallel plates 
[T]. Experimental advances in precision measurements 
of the Casimir force [5J [S] have revived interest in find- 
ing frameworks where one can compute these forces both 
numerically [H [5] and analytically. A particularly suc- 
cessful approach in applications to different geometries 
and material properties is based on scattering methods 
and techniques |6Vll2j. In this approach, the quantum- 
field-theoretic problem is reduced to that of finding the 
classical scattering matrix of each object. 

Another manifestation of fluctuations appears in the 
so-called dynamical Casimir effect: when objects are set 
in motion, they interact with the fluctuations of the back- 
ground vacuum in a time-dependent fashion which ex- 
cites photons and emits radiation. In fact, accelerating 
boundaries radiate energy and thus experience friction. 
An early example of this phenomenon was discussed by 
Moore for a one dimensional cavity [13]. A relativistic 
analysis of an accelerating mirror in 1-1-1 dimension in 
Ref. [14] employs techniques from conformal fleld the- 
ory. A perturbative study of the latter confirmed and 
generalized its results to higher dimensions [TS]. Among 
other methods, the fluctuation-dissipation theorem has 
been used to compute the frictional force on a moving 
sphere in free space |16j . a Hamiltonian formalism has 
been applied to the problem of photon production in 
cavities [T71 [TH] , and a (Euclidean) path-integral formu- 
lation is introduced to study the "vacuum" friction for 
a rough plate moving laterally (19J. We specially note 
that an input-output formalism relating the incoming 
and outgoing operators is used to compute, among other 
things, the frequency and angular spectrum of radiated 
photons [SO]- While a substantial literature is devoted 
to objects with perfect boundary conditions, dielectric 
and dispersive materials have also been studied in some 
cases [21]. In fact, dispersive objects exhibit similar ef- 
fects even when they move at a constant velocity. For 



example, two parallel plates moving laterally with re- 
spect to each other experience a (non-contact) frictional 
force [221 123] . Even a single object experiences friction 
if put in constant rotation 1 241 125) , a phenomenon most 
intimately related to superradiance flrst discovered by 
Zel'dovich [26]. (Translational motion of a single object 
is trivial due to the Lorentz symmetry.) The latter ex- 
amples, consisting of dispersive objects moving at a con- 
stant rate, are usually treated within the framework of 
the fluctuation-dissipation theorem or the closely related 
Rytov formalism [27) . 

Inspection of the literature on the dynamical Casimir 
effect leads to the following observations: There are 
a plethora of interesting — sometimes counter-intuitive — 
phenomena emerging from the motion of a body in an 
ambient quantum fleld [281 129) . These phenomena span 
a number of subfields in physics, and have been treated 
by a variety of different formalisms. Even the simplest ex- 
amples appear to require rather complex computations. 
Only recently experimental realizations — using a SQUID 
to mimic the moving boundary of a cavity (transmission 
line) [30l|3T] — have made precise measurements possible, 
raising the hope for an explosion of activity similar to 
the post-precision experiment era of static Casimir forces. 
This motivates reexamination of theoretical literature on 
the subject, aiming for a simple and unifying framework 
for analysis. 

In this work, we follow two goals. First, inspired by 
the success of the scattering-theory methods in (static) 
Casimir forces, we attempt at extending these techniques 
to dynamical Casimir problems. We find that the clas- 
sical scattering matrix is naturally incorporated into the 
formalism. However, dynamical configurations provide 
new channels where the incoming frequency jumps to dif- 
ferent values, hence the scattering matrix should be de- 
fined accordingly. Second, we aim for a universal frame- 
work which brings the diverse set of problems in dynam- 
ical Casimir under the same rubric. Most notably, we 
treat accelerating boundaries, modulated optical setups 
and moving dispersive objects — usually tackled with dif- 
ferent techniques, as explained above — on the same foot- 
ing. 

In this paper, computations are performed for a scalar 
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field theory. The generalization to electromagnetism is 
straightforward in principle, while practical computa- 
tions are more complicated in the latter. We find scalar 
field theory convenient to set the framework for more 
realistic applications. 

The paper is organized as follows: In Sec. |llj start- 
ing from a second-quantized formalism, we derive general 
formulas for the energy radiation due to the dynamical 
Casimir effect. In Sec. Ill we consider lossless objects 



undergoing non-uniform motion or optical modulation, 
and provide a variety of examples to showcase the power 
of the scattering approach. Specifically, we find that an 
(asymmetrical) spinning object slows down, and further 
contrast linear and angular motion. In Sec. |IV[ we con- 
sider dispersive objects moving at a constant rate. We 
also generalize to the case of multiple objects in relative 
motion where we study, among other things, an "atom" 
moving parallel to a dispersive surface. 



II. FORMALISM 

We start with input-output relations as described in 
Ref. [32]. The underlying formalism has been developed 
to quantize the electromagnetic field in a lossy or am- 
plifying medium |33H35| . (A similar method is also used 
to study the dynamical Casimir effect; see, for exam- 
ple, Ref. [20]. However the more general formalism in 
Ref. allows further extensions specially to dispersive 
objects.) Within this formalism the operators a'" and 
-out represent annihilation operators of the incoming or 
outgoing waves, respectively, in the vacuum (outside the 
object). These operators are then related by ^32 
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where b is the operator corresponding to the absorption 
within the object, and a and /3 are quantum numbers. In 
this equation, S is the object's scattering matrix while U 
describes its lossy character — the latter is related to the 
scattering matrix as shown later. This method treats 
field theory in the second quantized picture where quan- 
tum (annihilation or creation) operators are introduced. 
Equation ([l]) then relates quantum operators via the clas- 
sical scattering matrix. This proves to be useful in ap- 
plications to the dynamical Casimir effect. 

Equation ([T]) has its roots in the classical wave equa- 
tion. To see this, we first define "in" and "out" wave 
functions. The incoming wave can be expanded as 



$'"(x)=^C<"(x), 



(2) 



where c™ is the amplitude of the corresponding wave 
function The latter function should be normalized 
so that the number of incoming quanta per unit time is 
(negative) unity, 

jJd-E- [<i>r ^^n' - v$r = -s.^, (3) 



with the integral defined over a closed surface enclosing 
the object. The reason for this choice is that we shall as- 
sociate the wave function with a quantum operator which 
satisfies the canonical commutation relations, and the 
normalization should be defined consistently. Similarly 
the outgoing wave functions are normalized as 

^ ^ dS • [$°"** V$;^"* - V$°"** $;^"*] = (5„/3. (4) 

We also designate the solutions to the wave equation in- 
side the object as Now suppose that there are 
sources in two regions in space: at infinity where they 
generate the incoming wave, in Eq. (l2|); and within 
the object where they induce a field = z2a d-a^'^^ — 
the normalization of these functions does not affect the 
scattering matrix and thus is not discussed here. The 
incoming wave is scattered by the object while the ob- 
ject itself radiates due to the induced field. The resulting 
outgoing wave, = J2a cj^"**&a"*, is determined by 
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(5) 



where the first term is merely the scattering of the in- 
coming waves, and the second term captures the radi- 
ation of the object itself. The above analysis is based 
on a classical wave equation. Equation ([T]) extends the 
last equation to a relation between quantum operators, 
i.e. the complex-valued coefficients in Eq. ^ become 
quantum operators in Eq. ([T| through c™/°"' — > a'"/°"* 
and d ^ b (see also the discussion in Ref. i36j on the 
relation between wave mechanics and the classical limit). 
From this point on, we shall drop the hat symbol from 
quantum operators. 

For objects in motion, scattering can also change the 
frequency of the incoming wave. We make the depen- 
dence on frequency explicit while reserving a for other 
quantum numbers; the sum in Eq. ([T]) is replaced by 
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°° dw 
27r 



Most importantly, the scattering matrix may mix posi- 
tive and negative frequencies, in which case an outgoing 
operator of positive frequency is related to an incoming 
operator of negative frequency via Eq. ([T]). Note that an 
operator a^a with negative uj should be interpreted as a 
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creation operator; more precisely. Quia 
is related to a by time reversal. 

We assume that the environment is at a temperature 
Tonv while the object is at a (possibly different) temper- 
ature T. The distribution of the incoming modes (before 
scattering) is solely characterized by Tcnv, 

(aL"VL°a) = sgn(w)7i(a;, Tcnv) 6{uj - uj') S^p ■ (6) 



where n(a;, T) 



is the Bose-Einstein factor. 



cxp{hu]/kT)-l 

Note that this equation holds for both positive and neg- 
ative values of frequency; an operator a^a (ajj^,) defined 
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at a negative frequency, w, is interpreted as a creation 
(annihilation) operator of a positive-frequency mode. 

On the other hand, the occupation number for the op- 
erators b, locahzed on the object, is determined by T, 
the object's temperature. But we should keep in mind 
that the object could be moving, so the frequency defined 
from the point of view of a reference system (co)moving 
with the object is different from that of an observer in the 
vacuum, or the lab, frame. For a partial wave {uj,a) in 
the lab frame, we define Wq as the frequency according to 
the comoving reference frame. The occupation number 
is thei{3 



The choice of basis a is a matter of convenience, as in a 



Sgn(Wa)?l(tia, T) 6{uj - 

OOJ 



5^f3 ■ (7) 



From the distribution of the incoming and localized oper- 
ators, Eqs. (|6| and ([t]), we can evaluate the distribution 
for an outgoing mode (w', /3), 
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Sgn(w)7l(a;,ro„v) ^ |5'^'^,a;a|^ 



X! / -^^&^{l^a)n{,i^a,T)\Ua,'fi,uja\^ ■ (8) 



To find the flux of field quanta to the environment, one 
should compute the difference of outgoing and incoming 
flux. To study Eq. ([S]) in some detail, we consider two 
different situations. 

Accelerating objects — First we assume that the ob- 
ject is non-lossy so that the second term on the RHS of 
Eqs. ([l} and ([s]) is absent. For the sake of simplicity, we 
choose to work at zero temperaturej^ i.e. Tcnv = 0; the 
generalization to finite temperatures is straightforward. 
From Eq. ([s]), we find 
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We have used the fact that the Bose-Einstein distribution 
at T = is different from zero only for negative frequen- 
cies. Loosely speaking, this means that, in the vacuum 
state, all single-particle states with negative energy are 
occupied while those of positive energy are empty. The 
rate of energy radiation is obtained as an integral over the 
outgoing flux in Eq. ^ multiplied by the quanta energy 



(10) 



^ The change of basis from the frequency in the moving frame to 
that of the lab frame gives rise to the Jacobian. The partial 
derivative is positive on physical grounds. 

^ In the absence of loss, the object's temperature does not play a 
role. 



basis-independent notation Eq. ( 10 1 is cast as 



P = 



dw 
2 



(11) 



where S is the basis-free scattering matrix. These 
equations provide a simple and compact formulation 
which serve as the starting point for studying accelerat- 
ing boundaries and modulated optical devices in Sec. |III| 



Stationary motion — Next we consider objects in 
stationary, linear or rotational, motion. Although the 
objects are moving, the boundaries do not change their 
shape or orientation. One such example is two infinite 
plates moving parallel to their surface. Despite the mo- 
tion, the relative configuration of the two plates does not 
change in time. 

This type of dynamical problem is not explicitly time 
dependent, nevertheless the relative motion leads to dissi- 
pative effects. In other words, such systems respect time 
translation but break time-reversal symmetry, and thus 
allow for dissipation. Because of the stationary charac- 
ter of the setup, however, the scattering matrix § as well 
as the matrix U are diagonal in frequency, and we indi- 
cate this by a single frequency dependence as S'^a (w) and 

For a lossy object, the scattering matrix cannot be uni- 
tary as part of the incoming wave is lost inside the object. 
Interestingly, unitarity alone, sufficiently constrains the 
matrix U for our purposes [35]. There is a large body of 
literature on quantization in an absorbing (or amplifying) 
medium, covering a variety of approaches. The method 
of input-output relations [55HMj starts by formulating 
canonical commutation relations for the incoming and 
outgoing operators. 
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sgn{uj)5{uj - uj')6ai3- (12) 



We extend this method to moving systems by demanding 
that the operators b (localized on the object) satisfy the 
commutation relations in the rest frame of the object, 

[&c^Q,&L'^] = sgn(a)tt)^^5(cj- w')5„^, (13) 

with Cba defined above. This set of relations along with 
Eq. ([!]) lead to 

Oi Oi 

(14) 

In the comoving frame, the object is momentarily at rest, 
and thus the frequency according to this frame does not 



sgn(c.)(l-^|5^„(t.)|2)^^sgn( 



change, i.e. Wq = ujp if C/^q 0. Therefore, Eq. (14) can 
be recast as 

sgn(u;) {l-Y, \SUu:)\^) = sgn(^^)^ ^ \Up^{u:)\^ , 

a Q 

(15) 
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or in a matrix notation, 



sgn(c.) (I - S§t)^^ = sgn{Cjp)^{m^)p0 , 



(16) 



which constrains the matrix U in terms of the scattering 
matrix, §. In the hmit of static objects, one recovers the 
(basis-free) relation |32| 



I-§§^ =1LJU^. 



(17) 



This equation is interpreted by Beenakker as a 
"fluctuation-dissipation" relation with the LHS giving 
the dissipation due to the classical scattering from a lossy 
material, and the RHS accounting for field fluctuations 
due to spontaneous absorption or emission (in an ampli- 
fying medium) of field quanta. 

Equation (151 for a moving object can be inserted in 
Eq. ([S]) to obtain the flux due to the outgoing quanta. 
We are interested in the total flux, 
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The total radiation is obtained as the latter quantity mul- 
tiplied by hw integrated over frequency. Using Eqs. (|8|) 



and (151, one obtains the radiated energy per unit time 
as 

"P = / —hi^'^{n{[ba,T)~n{ljJ,Tcr,-v))[5l3a-\Spa{L0)\' 

(19) 



In the absence of motion [Cja = w), this equation cor- 
rectly reproduces the thermal radiation from an object 
out of equilibrium from the environment [321 137] . Inter- 
estingly, the moving object radiates energy even when 
the temperature is zero both in the object and the en- 
vironment. In this limit, the energy radiation takes the 
form 



= / ^^W^e(-(IiQ) [\Si3a{i^)\^ - Spa 

^ a,l3 



(20) 



where Q is the Heaviside step function. Therefore, spon- 
taneous emission takes place for a process whose fre- 
quency w is positive in the lab frame while, from the 
point of view of the moving observer, the corresponding 
frequency is negative. This mixing between negative 
and positive frequencies is at the heart of the dynamical 
Casimir effect 



In Sec. IV we employ Eq. ( 20 ) to flnd the spontaneous 



emission due to a rotating object. Furthermore, we study 
the conflguration of multiple objects in relative motion 
where we generalize the results presented in this section. 
In the process, we find that Eqs. (16) and (17) need to 
be modified for evanescent waves. 

In summary, Eqs. (10) and (20) express the energy 



moving object, respectively. In both cases, the radiated 
energy density is related to the off-diagonal part of the 
scattering matrix in Hence, they have a characteris- 
tic "Fermi-Golden-Rule" structure. In fact, to the lowest 
order, one can think of the off-diagonal S'-matrix as a 
potential due to the boundary condition or the object's 
material, akin to the Fermi Golden Rule. 

In the following sections, we provide a variety of exam- 
ples where we discuss applications of the general formulas 
presented above. 



III. LOSSLESS ACCELERATING OBJECTS 

In this section, we consider lossless objects with differ- 
ent shapes in various dimensions undergoing rotational 
or translational motion or oscillation. In the process, 
we reproduce some existing results in the literature, and 
also present many novel applications. Equation ( 10 ) is 



the central formula according to which we compute and 
discuss these results. 



A. A Dirichlet point in 1 + ld 

The prototype of dynamical Casimir phenomena is the 
motion of a point-like mirror in one dimension |141 115j . 
For simplicity, we assume that the ambient vacuum con- 
sists of a scalar field, $, subject to Dirichlet boundary 
conditions on the mirror 



$(t,g(t)) = 0, 



(21) 



where q{t) is the trajectory of the mirror in time. We use 
a perturbative scheme |15) where we expand Eq. (21 ) for 
small q{t) to obtain 



^t,o) + q{t)dMt,o) 



0. 



(22) 



The scattering solution can be formally expanded in pow- 
ers of q, 

$ = $0 + *i + • • • • (23) 
The boundary condition, to the first order, takes the form 
$i(t,0) = -g(t)9,$o(i,0). (24) 
The incoming and outgoing modes are defined as 



^in/out 



- exp[—iijj{t ± z/c)]. 



(25) 



Note that the normalization is chosen in accordance with 
Eqs. ([3]) and In the zeroth order, i.e. for a static 
mirror, we have 



(26) 



radiation for an accelerating lossless body and a lossy 



We can compute $1 by solving the free fleld equation 
(□ $ = 0) in the vacuum subject to its time-dependent 
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value at the origin (z = 0) as given by Eq. ( 24 ) . We leave 



the details to Appendix |A 1[ the scattering matrix (from 
either side of the point) is obtained as 



2iq{n) 



(27) 



where q{i^) is the Fourier transform of q{t). (Note that, 
here and in the following, we only write the off-diagonal 
correction to the scattering matrix.) One can then com- 
pute the radiation according to Eq. (10 1. For a Fourier 



mode f2, the integral in the latter equation contributes 
in the window of < — cj < il. Putting all the pieces 
together, and multiplying by a factor of two accounting 
for the scattering from both sides, we find 
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(28) 



In the last line, the radiation is expressed as an integral 
over time. From Eq. (28), one can infer the dissipative 
component of the force 



fit) 



(29) 



in complete agreement with Ref. [15j 



B. Modulated reflectivity in 1 + ld 

For moving bodies, the dynamical Casimir radiation 
is difficult to detect experimentally since it requires the 
objects to move at very high frequencies. An alternative 
approach is suggested by modulating optical properties 
of a resonant cavity [351 131] • In fact, any linear time- 
dependent process can lead to similar dynamical Casimir 
effects. Modulated reflectivity, for example, generates 
photons and gives rise to radiation [28l |40] . The latter 
can be studied within the same framework that we devel- 
oped for non-lossy objects]^ In this section, we consider 
a point particle in one spatial dimension, but, unlike the 
model in the previous (sub)section, a linear coupling of 
(time-dependent) strength e is introduced at the position 
of the particle. The field equation then reads 



1 



d^A^{t,z) + e{t)S{z)^{t,0) ^0. (30) 



We recover a perfectly reflecting object for e — ?> 00. An 
imperfect mirror undergoing arbitrary motion is studied 



^ In treating the dynamical Casimir effect in the absence of loss, 
we did not assume that the objects are actually moving. 



in Ref. [H]. Note that in our model the particle is at 
rest at the origin while the coupling is modulated. The 
5-matrix can be computed by techniques similar to quan- 
tum mechanical scattering in a one-dimensional delta po- 
tential. For simplicity, we take e(t) — eo+eo cos(r2i) with 
Co ^ esi- We note that there are new scattering chan- 
nels with incoming waves from one side transmitted to 
the other side of the object. A scattering ansatz with 
incoming waves from the RHS is given by 



f (j)-Rin ^_ out 



z > 0, 
z < 0, 



(31) 



where summation is made over both signs, and the wave- 
functions denoted by R{L) are defined on the right (left) 
side. Obviously, the two sets of definitions are related 
by reversing the sign of the coordinate z. In the above 
ansatz, we have exploited the smallness of the oscillatory 
part of e by truncating the sum at the lowest harmonics. 
One can obtain the scattering amplitudes by matching 
the functions on the two sides of the mirror while setting 
the difference in their first derivative to e$(t, 0). We find 



r± = ti 



ien ^|w(a; ± fi)]/ c 



(eo - 2iw/c)(eo - 2i(w ± fi)/c) 



(32) 



This equation can be further simplified by assuming eo ^ 
ri/c. In this limit, the energy radiation (per unit time) 
is, according to Eq. (10), 



(33) 



For a general e{t) slowly varying around a mean value of 
eo, the total energy radiation is given by 



P = 



(34) 



C. A Dirichlet line in 2+ld 

Now consider a line extended along the x axis in two 
spatial dimensions. In this geometry, the incoming and 
outgoing waves are described by 



n./out 1 



: exp(— icji -|- ikxX =p ikj_z) , 



(35) 



where is the wavevector along the line while kj_ is 
the perpendicular component, k±{uj,kx) — y^uj^/c^ — k"^. 
We assume that the line undergoes a rigid but time- 
dependent motion, q{t), normal to the x axis. One then 
finds 



Suj+nk^,ujk^ = -2iq{p.)^k^{i^,kx) kj_{u 



(36) 
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see Appendix |A 1| The scattering matrix is diagonal in 
kx due to translational symmetry along the cc-axis. Note 



that Eq. ( 36 1 is computed only for propagating modes- 



evanescent waves fall off rapidly with the distance from 
the surface and do not contribute to radiation at infin- 



ity. The sum over all partial waves in Eq. ( 10 ) becomes 
/ "^2-^°' with L being the extent of the line; the integra- 
tion is over propagating waves only, i.e. c\kx\ < \uj\ and 
c\kx\ < u + Q. Finally, an integration over the frequency, 
w, gives 



P = 



hL 

128c3 



2^ 



g(r!)|2r!5 



(37) 



In order to write this equation in the time domain, we 
extend the integral to (— oo, oo) and recast the integrand 
as SI |g(r2)p Imx(f^). The function x is the (normal- 
ized) response function whose imaginary part is pro- 
portional to the energy dissipation to the environment, 
Imx(S^) — |ri^sgn(r2). Because of causality, the full re- 
sponse function can be obtained via Kramers-Kronig re- 
lations. In the time domain, we find this function as 



24 1 



with P being the principal part. In this context, the 
response function relates the force to the object's dis- 
placement. So the force acting on the object at time t is 
given by 



fit) 



3hL 



dt'P 



1 



{t-t'f 



q{t'). 



(39) 



The force is manifestly causal, i.e. it depends on q at ear- 
lier times. However, Eq. (39) is possibly divergent near 



the upper bound of the integral unless a short-time cutoff 
is introduced to replace this bound by t—r. This does not 
affect the dissipative component of the force but regular- 
izes the inertial force which sensitively depends on the 
large-frequency physics. In fact, in deriving Eq. (39), we 



have used the small-frequency behavior of the response 
function, so this equation should be valid only for large 
times. Interestingly, the force does not vanish even when 
the object no longer moves. We find that the force, long 
after the object comes to a full stop, falls as 



fit) 



3hL 



dt'qit')^^m, (40) 



where the integral is over the displacement of the object 
for the duration of the motion (which is assumed to be 
much smaller than t) , and q{0) is the integrated displace- 
ment. 



Figure 1. A segment in a waveguide. The arrows indicate the 
direction along which the segment oscillates. Below a certain 
frequency, cjmin = 2cjo, the motion is frictionless. 



(a waveguide) in two spatial dimensions; see Fig. [T] 
Dirichlet boundary conditions are assumed on all sur- 
faces. The segment undergoes a rigid motion q{t) parallel 
to the infinite lines. The only difference compared to the 
previous (sub)section is that the modes along the x-axis 
are quantized. Therefore the integral over k^ is replaced 
by a sum over n where kn = The radiation takes the 
form (cf. Eq. (l37|) 



P = 



128c2 



<m 

2^ 



q{n)\^^^9i^L/c), (41) 



where the function g is defined as 



(38) g(z.) = 



512 



E 



dz(l-z)W(2 



-)((1-^)^- 



-)• (42) 



For large L, g{v) — > and we recover the results of 
the previous (sub)section. However, this function van- 
ishes below 1/ = 27r; a low- frequency motion does not 
dissipate energy since propagating waves inside a waveg- 
uide have no support in the range (— wqjWo) with 
being the lowest eigenmode of the waveguide, hence 
<^min = 2a;o = I'KcjL. Close to this frequency, the func- 
tion g vanishes quadratically. 



9ii^) 



V > 2tt. 



(43) 



Similar to the previous (sub)section, one can use 
Kramers-Kronig relations to obtain the response func- 
tion. Specifically, we are interested in the long-time limit 
after the object comes to a full stop. The dependence on 
large t can be inferred from the short-frequency response 
(Eq. (HI)) as 



9(2c^o)), 



(44) 



with ujQ = TTc/L as defined before. Note that the force 
now falls as 1/t^ while its amplitude undergoes periodic 
oscillations at frequency 2ajo, twice the lowest natural 
frequency of the waveguide. 



D. A Dirichlet segment in a waveguide in 2+ld 



E. A Dirichlet plate in 3+ld 



Next we consider a finite segment of size L along 
the X-axis confined between two infinite Dirichlet lines 



In this section, we consider a (two-dimensional) plate 
in three dimensions subject to the same (Dirichlet) 
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boundary conditions. Again we assume that the plate 
undergoes a rigid motion q{t) normal to its surface, and 
find the scattering matrix as 



-2i q{n) ^ kj_{uj,k\\) kj_{uj + n,^^) ; 



see Appendix |A 1[ In order to compute the radiation, 
one should integrate over propagating modes only (both 
for incoming and outgoing waves), i.e. c|k||| < |a;| and 
c|k|| I < oj + n. Equation ( 10 1 then gives 



1807r2c4 



3607r2c4 



2^ 



dtq- 



(46) 



where is the area of the plate. The resulting (dissipa- 
tive component of the) force is then 



m = - 



3607r2c4 



(5) 



(47) 



again in agreement with Ref. |15j . 

We also note the difference between odd and even di- 
mensions. In 2 dimensions, the force displays long-time 
tails, while in 1 and 3 dimensions it is an (almost) in- 
stantaneous function of the displacement. 



F. A Dirichlet corrugated plate in 3+ld 

We can generalize the results in the previous 
(sub)section by considering a corrugated plate. For cor- 
rugations of wavevector K, the scattering matrix is given 
by (cf. Eq. (|45|) 

5'tj+ak||+K,ijk|| = 



-2iq{n,K)^kj_{uj,k\\)kj_{uj + n,k\\ +K). (48) 

where q, the displacement from the x — y plane, is a 
function of both w and K; see Appendix |A 1[ The con- 
dition for propagating waves is modified as c|k|j| < \uj\ 
and c|k|| -I- K| < uj + il. The radiation formula should 
be modified accordingly to ensure that only the prop- 
agating modes are integrated. While this integral can 
be computed explicitly, we use a trick as follows: For 
f2 > c|K|, we Lorentz-transform to a frame in which 
the scattering matrix is diagonal in the wavevector ky; 

the velocity of this frame is v = The scatter- 

ing matrix, the integral measure, as well as the con- 
dition for propagating waves are invariant under such 
a transformation. But the frequency {huj' in Eq. (10)) 



picks up a factor of 7(v) — — v'^/c^, while the 

lower bound of the integral over uj changes from — O to 
— 7(v)(fi — V • K) = — il/7(v) which, through compar- 
ison with the first line of Eq. (|46|), contributes a factor 



of 1/7^. Hence, the radiated energy density in and K 
becomes 

Pin^K) ^ ^^^\q{n,K)\'n{n' ^ c'K')'/^ (49) 



Note that the difference of a 



consistent with Ref. [ ^ 

factor of two in comparison with Eq. ( 46 ) is in harmony 
with the setup in Ref. ^Vd^ where the plate occupies a 
half-space. Similar results can be obtained for Neumann 
boundary conditions |42j . 



G. A Dirichlet sphere in 3+ld 

In this section, we consider a sphere subject to Dirich- 
let boundary conditions linearly oscillating in three di- 
mensions. The oscillation amplitude, q{t), is small com- 
pared to the radius of the sphere, R. We choose the z 
axis parallel to the motion and passing through the cen- 
ter of the sphere. The incoming and outgoing waves for 
a spherical geometry are defined as 



in/out 



oj/m 



where ft-^'^'' are spherical Hankel functions, and Yim is 
the usual spherical harmonic function. Due to azimuthal 
symmetry, the scattering matrix is diagonal in the index 
m but possibly mixes different Is. The scattering from 
an oscillating sphere can be computed by using Green's 
theorem; see Appendix | A 2| for more details. We find the 
(off-diagonal) scattering matrix as 




''^^ ['^)yUO,^). (50) 



LO~\-Ql'm., Lolm 



2iq{n) 



du'r 



X y/\{uj + il)uj\Fi 



ujR 



{Lo + n)R 



(51) 



where du'm as defined in Appendix |A 2| is nonzero only 
for /' = i ± 1, and the function F is defined as 



1 



xh^\x) 



(52) 



Using Eq. (10), the radiated energy density is given by 

8h\q{n)\'^ 



p{n) 



3c2 

OO 



do; 
2^ 



1=0 



(53) 



We consider two different limits: 

a) Q.R/c ^ 1. For a slowly oscillating object, we need 
only consider the lowest partial wave, i.e. the I = term 
in Eq. ( |53p . One then finds the radiated energy density 
in frequency as 



p{n) 



hVL^R^ 
307rc4 



(54) 
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b) ^R/c ^1. In this case, one should include all 
partial waves up to ^max ~ flR/c ^ 1. Below we closely 
follow the line of argument in Ref. [43,. For large I, we 
have 



\Fiix)\^il-lVxy/\ 



1 < / < X. 



(55) 



The sum over all partial waves can then be recast as an 
integral over / yielding (through the change of variables 
a = l/{nR/c) and x = \uj\/n) 



4nR^n'^\q{n)\ 



1/2 

daa 



dx (1 — x)x'^ 



1 - 



1/2 



1 - 



(1 



1/2 



hfi^R^ 

2707rc4 



\m)\' 



(56) 



This equation reproduces the contribution of the TE 
modes to the electromagnetic version of aperfectly re- 
flecting sphere; see Eq. (4.20) in Ref. [ISjF] Indeed one 
finds a similar correspondence for a perfectly reflecting 
plate, that is the radiation due to the TE modes is equal 
to that of the Dirichlet plate [20] . 

Finally we note that the radiation due to the oscillatory 
motion can be computed for a variety of other geometries 
such as cylinders, ellipsoids, etc. 



H. A spinning object in 2+ld 

Heretofore, we studied examples where the object is 
accelerated by an external force. In this section, we con- 
sider an object rotating at a constant angular velocity 
ri. We further assume that the object is not rotationally 
symmetric, as illustrated in Fig. [2j As usual, we assume 




Figure 2. An (asymmetrical) spinning object. The object 
slows down as it emits "photons". 

a scalar field subject to Dirichlet boundary conditions on 
the object's surface, and limit ourselves to 2-f 1 dimen- 
sions. 

Since the orientation changes with time, waves imping- 
ing on the object are partially scattered at a shifted fre- 
quency determined by Vl. The scattering matrix is more 
conveniently computed by going to the object's reference 
frame. 



There are however mixed terms between the two polarizations 
which vanish in this limit, see Ref. |43|. 



I 

We first start with the field equations in the laboratory 
(static) frame. The wave equation, □$ = 0, in polar 
coordinates is 



\dt - -drrdr 
r 



$(t,r,0) = O. (57) 



The incoming and outgoing waves are defined as 



^in/out^. 



(58) 



where H is the Hankel function. (We have dropped an ir- 
relevant constant in the definition of the these functions.) 
The rotating frame is described by 

t' ^t, r' = r, 4)' = (t>- fit. (59) 

The field equation in the latter frame takes the form 



1 



1 



^{dt - fld^'Y drvdr 



1 



$'(<,r,0')-O. 



(60) 

Note that $'(i,r, 0') = ^{t,r,<j)). Specifically, in the 
new coordinate system, the functions ^i^rn as defined in 
Eq. (58) become 



K-nmmit,r,cf>') = e-(-f^™)*e'™^'H(i.2) {^) . (61) 

The rotating frame is more convenient to write the scat- 
tering ansatz as the object does not move in this frame, 
and thus the time dependence drops out as a phase fac- 
tor. In the latter frame, the boundary conditions take 
the form 



c 



c ^^^'<^>' H-(i) f - ^{m ~ m'))r 



0, (62) 



where E denotes the boundary. One can see that this 
equation indeed satisfies Eq. (|60]) once the time depen- 
dence, g-^l^-f^™)*^ is restored. The scattering matrix 
sends the frequency a; to a; — 17 (m — m!) from the point 
of view of an observer in the lab frame. To obtain an 
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analytical expression for the scattering matrix, we con- 
sider the non-relativistic limit where the object's (lin- 
ear) velocity is small compared to c. It then suffices 
to compute the scattering matrix for the lowest partial 
waves. As a specific example we consider an ellipse close 
to a circle of radius R with 5 being the difference of 
the two semiaxes, i.e., in polar coordinates, defined as 
r(0) = R + cos(2(/))(S/2. The scattering matrix to the 
lowest order in 5 is then obtained as 



iiriuj + 2n)^RS 
8c2 log(|w|i?/c) ' 

i'KLo'^RS 
8c2 \og{{ij + 2n)R/c) 



(63) 



The energy radiation per unit time, according to Eq. ( f 1, 
is 



nhR^S^n^ 
I0c4log(f2i?/c)2 ■ 



(64) 



Therefore, an (asymmetrical) object which is spinning, 
even at a constant rate, slows down due to quantum dis- 
sipation. Note that the torque (due to the back-reaction) 
is simply the radiation rate divided by the frequency il. 



I. A Dirichlet disk in 2+1 dimensions: linear vs 
angular motion 

Now consider a circular disk of radius R subject to 
Dirichlet boundary conditions in two spatial dimensions. 
Below we contrast two different types of motion, see 
Fig. [3] First we consider a linear oscillation, q{t) = 




Figure 3. Linear vs angular motion; the radiated energy is 
comparable in the two cases. 

S COS fit (with 6 being the amplitude of the oscillation) 
along the x-axis. The scattering from an oscillating disk 
can be obtained by using Green's theorem; see Appendix 
A 3 for more details. We obtain the scattering matrix as 



9 - ^ M 

ttR 



where the function M is 



ujR 



M, 



m± 1 



(cj + n)R 



(65) 



At low frequency, ^R/c ^ I, the scattering matrix for 
the lowest partial waves is obtained as 



Ac \og{ujR/c) 



(67) 



The energy radiation (per unit time) can be computed 
from Eq. (10); for a linearly oscillating disk, we find 



p. 



t 



64c2 \og{nR/cy 



(68) 



Next we consider the same disk undergoing orbital mo- 
tion. The latter can be thought as a rotation around a 
point off the center of the disk. Specifically, we assume 
that the disk's center undergoes a trajectory (r, (j)) — 
(S, fit) in polar coordinates, i.e. the center is at a fixed 
radius 6 while orbiting around the origin at frequency Q. 
The scattering process can be examined by techniques 
similar to the previous (sub) section where we considered 
a spinning object. We must assume ilR/c <C 1 as well as 
6 ^ R — similar to the linear motion. The latter allows 
us to compute the scattering matrix only for the lowest 
partial waves, as 



'S'w-i-ni.Ljo 

•Std-l-ao,!.!; -1 



Z7r(cj -|- il)S 
~2clog(|w|i?/c) ' 

inujS 
2clog{{uj + n)R/c) 



(69) 



Note the similarity between Eqs. (671 and (69). The 



radiated energy by a revolving disk can then be computed 
according to Eq. (10), as 



24c2 log(f7i?/c)2 ■ 



(70) 



It is interesting to compare Eqs. (68) and (70 1 where 



linear and angular motion have been considered, respec- 
tively. We notice that, apart form a proportionality con- 
stant, the analytical form of the dissipated energy is iden- 
tical in the two cases. 

The above results can be compared with the previous 
(sub)section. In the former case, the setup is symmet- 
ric under inversion with respect to the origin while the 
latter is not. Thus, in the lowest order, scattering of 
waves changes the angular momentum by two units for 
a spinning ellipse and by one unit for a disk in (linear 
or circular) motion; cf. Eqs. (63), (67) and ([69|. Conse- 



quently, the energy radiation in the former case, Eq. (64 ), 



is suppressed by two orders of magnitude in comparison 
with Eqs. (lesl) and ([70| for a disk. 



IV. 



STATIONARY MOTION OF LOSSY 
OBJECTS 



(66) 



In this section, we consider lossy objects but limit 
ourselves to constant (linear or angular) speed. Among 
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other results, we reproduce existing formulas in the liter- 
ature with significantly less labor. Our starting point is 
Eq. (pO| in application to rotating objects. We then gen- 



eralize to the case of multiple objects in relative motion. 



A. Rotating object 

Let us consider a solid of revolution spinning around 
its axis of symmetry at frequency fl. We choose polar 
coordinates (r, 0, z) where the z direction is along the 
axis of symmetry. The latter coordinates describe the 
lab frame in which the object is rotating. The rotating 
(or comoving) frame is defined by the coordinate trans- 
formation 



t' = t. 



r, (f)' 



nt, 



(71) 



A partial wave defined by frequency uj and azimuthal 
index m in the lab frame is characterized by the frequency 



= UJ — ilm 



(72) 



from the perspective of the rotating frame; see Ref. [5S] 
for a detailed discussion. In harmony with the discussion 
in Sec. |llj the object spontaneously emits energy when 
UJ > and ujm < 0, i.e. in the frequency window 



< UJ < ilm, 



(73) 



or the so-called superradiating regime first introduced by 
Zel'dovich |26| . Therefore, the energy radiation per unit 
time from a rotating object to the environment at zero 
temperature is given by 



V 



duj 
2ti 



huj Tr 



e{ni^~-uj) (§(w)§(a;)^ - I) 



(74) 

where Iz is the z-component of the angular momentum 
operator in units of h. Note that the scattering matrix is 
diagonal in frequency, uj, since the object is undergoing a 
stationary motion with its shape and orientation fixed in 
time. Equation ( 74 1 indeed gives the spontaneous emis- 
sion by a rotating object consistent with Ref. [35] where 
the Rytov formalism is used through an involved anal- 
ysis. For a small object, only the lowest partial waves 
contribute to the radiation, and we recover the results of 
Ref. [21] . Note that in deriving Eq. ( 74 ) we did not use 



any approximations regarding the velocity of the rotating 
object. 

For the sake of generalization to multiple objects, we 
point out that Eq. ([T]) can be interpreted in a simple way 
to arrive at the same results. According to this equa- 
tion, the probability amplitude for spontaneous emission 
is given by 



(75) 



where we have suppressed all quantum numbers other 
than UJ and m, and used the fact that the amplitude 



is diagonal in m due to the rotational symmetry of the 
object. The rate of this process is 



1, 



(76) 



where, in the last equality, we have used Eq. (151 in the 
superradiating regime, i.e. for < w < flm. Note that, 
in this regime, [^^((u;)! > 1, hence superradiance. The 
integral of Af multiplied by Huj (over superradiating fre- 
quencies) reproduces the energy radiation as given by 
Eq. (iTl. 



B. Moving plates 

A system comprising two lossy parallel plates under- 
going relative lateral motion is the canonical example of 
non-contact friction. In the following, we sketc h a simple 
derivation of this friction based on Eqs. ^ and ([TtIP] We 
first note that Eq. (17) can be interpreted via a classical 



argument. One can denote the RHS of this equation as 
the rate of "photon" absorption in a dispersive medium. 
Current conservation implies that the latter should be 
equal to the influx of the field quanta outside the body. 
Let us consider a classical wave scattered from the object 
as 



/3 



(77) 



where we have suppressed the frequency, w. The 
current density going into the body is given by 
=1 [$* v$ — V^* $]. Since $Q, and $^ are properly nor- 
malized (see Eqs. ([s]) and ^ and the explanation there- 
after), the total influx of fleld quanta is 



I 

2 / 



(78) 



But this is exactly the LHS of Eq. (17) 



To study moving plates we need to extend Eq. (17) to 



evanescent waves which arise in non-compact geometries 
(plates, cylinders, etc.), but are absent for a compact ge- 
ometry. Since such waves are not propagating, "incom- 
ing" and "outgoing" wave functions lose their straight- 
forward interpretation. In other words, they do not carry 
currents 



(j)Out/in* y^out/in _ ^^out/in* ([)Out/in 



0, 

(79) 



but satisfy a different relation (after proper normaliza- 
tion) 



dS- 



(80) 



^ For a uniform translational motion, one can Lorentz-transform 
Eq. | |17^ to the moving frame. 
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Therefore an incoming wave scattered from the object, 
$ = $Q + J2i3 Spa^"^'^^, carries an influx of "photons" 
given by 



(81) 



of current dictates 



where Eqs. ( 79 ) and ( 80 1 are used. Then the conservation 



(82) 



Equations (17 1 and 182^ define the matrix U in terms of 



the scattering (or reflection) matrix for propagating and 
evanescent waves, respectively. 

We can now compute the friction between two "dielec- 
tric" plates moving in parallel. We assume that the first 
plate is at rest while the other plate, separated from the 
first by a distance d, moves at a constant velocity v along 
the X axis. Because of translational symmetry, all matri- 
ces are diagonal in the frequency w and the wavevector 
k|| parallel to the surface. Here, Eq. ([Tj) finds a two-fold 
application. On one hand, it allows for spontaneous emis- 
sion from an object, while on the other hand, it describes 
the reflection and absorption of waves by a second object. 

1) Spontaneous emission; the source fluctuations in the 
first plate give rise to outgoing wave fluctuations. The 
amplitude of the spontaneous emission is given by 



(83) 



where the dependence of the matrix [/ on cj and k|| is 
implicit. Note that incoming waves do not contribute to 
spontaneous emission. 

2) Reflection: These outgoing waves propagate to 
the second plate and get a factor of e**^^'^ with k± — 
^(jj'^/c^ — ky ; a phase factor for propagating waves while 

exponentially decaying for evanescent waves. There they 
are partly reflected and partly absorbed by the second 
plate. The amplitude for "photons" spontaneously emit- 
ted by the first plate and then absorbed by the second 



A2^i = e 



ik\ d 



U2U1. 



(84) 



Equivalently, the rate of the latter process is given by 



ni 



ik\ d I 



\U2nUi\'ni, 



(85) 



where ni — n{uj,Ti) is the Bose-Einstein occupation 
number defined at temperature Ti. The superscript 1*** 
indicates that Eq. ( [85| is computed within the first re- 
flection. One can similarly compute Mi^2, the current 
from the second to the first plate. But in the latter case, 
n2(w,k||) is centered at uj — vkx, i.e. n2 = n{u; — vkx,T2), 
because the thermal fluctuations are defined with respect 
to the comoving framej^ The total flux from the first to 



For relativistic velocities, one should also include the Lorentz 



factor 7(1)) = l/\/l — jc? 



the second plate, within the first reflection, is then 



-'^2^1 



-Mll2^V''A\U2?\ViWn^-n2). (86) 

One can easily sum the contributions from multiple re- 
flections, 



n=0 



(87) 



1 - e^'i'^'^RiR2 ' 

where i?i and R2 are the reflection matrices. Note that 
the n-th term in the last equation is the amplitude for a 
"photon" spontaneously emitted by the flrst plate (C/i), 
reflected n times from the two plates {{e'^^^-^'^ R1R2)'"') 
before finally getting absorbed by the second plate (C/2)- 
The amplitude A\°^2 obtained similarly. The total rate 
then becomes 

A/2^1 - A/1^2 = -h _9,t., rfn n 19 ("1 ~ "-2)- 
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Also note that, from Eqs. (17) and (|82|, we have 

1^ |2 _ -L - i^tji , for propagating waves, 
" for evanescent waves. 




(89) 



/ 



Finally friction is the rate of (lateral-)momentum transfer 
integrated over all partial waves, 

duj f dk|| ^, |e*'=^'*ht/2p|f/i|2 

2^jj2W^^^^ \l^e^^^^<^R,R2V ~ 

(90) 

We should emphasize that the reflection matrix for the 
second plate should be computed in its rest-frame, and 
then transformed to the lab frame according to Lorentz 
transformations. The last equation is the analog of the 
results in Refs. (22l [23] for the scalar field. 

To be more concrete, we consider a scalar model that 
is described by a free field equation in empty space while 
inside the object a "dielectric" (or, a response) function 
e is assumed which characterizes the object's dispersive 
properties. The field equation for this model reads 

.2 , r72^ 



(e(a;,x)a;2-HV2)$(a;,x) =0, 



(91) 



with e being 1 in the vacuum, and a frequency-dependent 
constant inside the object. 

For a semi-infinite plate, the reflection matrix R is 
given by 



Ru'k.. — 



(.LiP- jf? — kn — .IlJ^Ic? — k' 



ea;2/c2 - kjj + . hJ^ j - k' 



(92) 



This is easily obtained by solving the field equations in- 
side and outside the plate and demanding the continuity 
of the field and its flrst derivative along the boundary. In 
a moving frame, the frequency and the wavevector should 
be properly Lorentz-transformcd. 

These reflection matrices can then be inserted in 



Eq. ( 90 ) to compute the frictional force. 
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C. An atom moving parallel to a plate 

In this section, we consider a small spherical object, 
an atom, moving parallel to a plate. In the non-retarded 
limit, an electrostatic computation is done in Ref. [53] for 
a similar setup. For our purposes, it is more convenient 
to consider the rest frame of the "atom" in which the 
plate moves laterally. This is another example of station- 
ary motion where the geometrical configuration does not 
change even though the objects are undergoing relative 
motion. We assume a small spherical object of radius 
a (much smaller than the separation distance d), such 
that the first-reflection approximation suffices. Similar to 
the previous (sub)section, we first consider spontaneous 
emission by each object. The plate (denoted by sub- 
index 2) emits "photons" of frequency uj and wavevector 
kii with a probability amplitude 



^2(w,k||) = [/2(cJ,k||). 



(93) 



Then, these waves propagate to and reflect from the 
"atom." Planar waves pick up a factor e'^'^^'^ upon travel- 
ing a distance d. To find the scattering off of the spherical 
object, we must change to a basis of spherical waves. A 
planar wave can be expressed as a superposition of spher- 
ical waves as 



47r 



z'j,(c.r/c)y,„(x)y;„(k) 



(94) 



where x and k are the unit vectors parallel to the vectors 
X and k, respectively. A planar wave, ^wky , defined with 
respect to a reference point on the plate's surface below 
the sphere's center is related to spherical waves, ^uhm 
centered around the "atom" as 



<j)OUt 

^wkii 



^ 27rz'e»^-^'^l^;,(k) 

Im 



Kt^) ; (95) 



see Appendix \K\ for the definition of planar and spherical 
waves. Then the amplitude for "photons" spontaneously 
emitted by the plate and then absorbed by the sphere is 



^2 



C/i,,„(^)t/2(c^,k,|), (96) 



where Ui characterizes the loss due to the "atom." Sim- 
ilarly, we can compute the amplitude .42<-i for the in- 
verse process where the spontaneous emission due to the 
"atom" is absorbed by the plate. One can then obtain 
the rate of energy or momentum transfer between the 
objects. An analysis similar to the previous (sub)section 
gives the force within the first reflection, 



/ 



^ik\ d I 



|11^(k)P|?7i;„(^)P|£/2(^,k||)|- 

|fc_L|a;/47r2c 



(97) 



where ni(a;) = n{oj,Ti) and 7x2(0;, k) = n(w — vkx,T2) 
are the Bose-Einstein factors for the "atom" at temper- 
ature Ti and the plate at temperature T2, respectively. 
Note that we have only considered the flrst reflection as 
the "atom" is small compared to the separation distance. 
The matrix U2 is given as in Eq. ( 89 1 while, for the spher- 



ical object, there is no evanescent wave and thus the ma- 
trix Ui is constrained by 



|C^ii,„(a;)|^ 



1 



\SUuj)\^ 



(98) 



with Sim{<^) being the scattering matrix of the "atom." 
This equation indicates that a frictional force (or en- 
ergy transfer) arises only if the "atom" is lossy, i.e. 
[^(^(a;)! < 1. For the scalar model introduced in the 
previous (sub)section, the sphere's scattering matrix is 



hf'\uja)daji{n{uj)u)a) — dah['^\u}a) ji{n{uj)uja) 



hf\oja)daji{n{ij:)uja) - dah\'^' (uia) ji{n{u))uja) 

(99) 



(1)/ 



where a is the sphere's radius, and n(a;) = a/ £5(0;) with 
es being the "dielectric" function of the spherical object. 
To the lowest order in a/d, we shall limit ourselves to the 
low-frequency scattering of the partial wave I — m — oj^ 
Within this approximation, the friction at zero tempera- 
ture takes the form 



/ 



dkii 



fex>0 



duj — ^ 







(100) 



The reflection matrix R can be obtained from Eq. (92) 
via Lorentz transformation. Similarly, one can consider 
the frictional force between a rotating sphere and a sta- 
tionary plate [IS]. With our scalar model, the scattering 
matrix for a rotating sphere is obtained from Eq. ( 99 ) 



by changing the argument of the Bessel functions to 
n{u!„i)uj,na where = uj — Om. Having the scatter- 
ing matrices of a moving plate and a rotating sphere, 
one can compute the friction when both objects are set 
in motion. 



V. SUMMARY AND OUTLOOK 

In this work, we have developed a unifled scattering 
approach to the dynamical Casimir effect. We have ob- 
tained general formulas for the radiation from moving ob- 
jects for accelerating boundaries and modulated optical 
devices without loss, as well as lossy bodies in uniform 
motion. We provided and studied numerous examples, 
many of which are novel, to better illustrate the techni- 
cal power and conceptual elegance of the scattering ap- 
proach. 



^ One should note that in the case of electrodynamics there are no 
monopole fluctuations and thus the leading contribution to the 
friction comes from 1 = 1 |44| . 
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Extensions to more realistic boundary conditions and 
field theories should be of both theoretical and practi- 
cal interest. A superfluid liquid, for example, presents 
a natural framework to study the effects of motion in a 
quantum vacuum. Indeed, similar effects such as radia- 
tion and spontaneous emission have been discussed for a 
superfluid [15] . 

Quantum electrodynamic requires a further treatment 
due to its vector nature. An obvious complication arises 
since electromagnetic waves are polarized, enlarging the 
scattering matrix, and complicating practical computa- 
tions. Realistic materials also include loss. In Sec. |III[ 
we considered accelerating objects with perfect boundary 
conditions. Extensions to lossy objects require a general- 
ization of the formalism presented in Sec. |TTj An expres- 
sion solely in terms of the scattering matrix is desired in 
the latter case when object's acceleration and dispersive 
properties interplay in a rather complex fashion. One can 
anticipate that the computation of the scattering matrix 
will be more complicated in this case. It is also worth- 
while to consider configurations of multiple objects in 
arbitrary motion. We have partially tackled this prob- 



lem in the context of stationary motion in Sec. IV, while 
extensions to accelerating objects will present new chal- 
lenges and provide further insights. Specifically, one can 
ask how the (inertial as well as dissipative) forces be- 
tween two objects change as the result of their motion or 
acceleration. 

The formulation of dynamical Casimir effect in terms 
of the scattering matrix should also provide an efficient 
prescription for numerical computations. The scattering 
matrix is purely a classical quantity, and presumably can 
be numerically computed with high precision. This is 
particularly important if the motion cannot be treated 
perturbatively — when the speed, the amplitude of oscil- 
lations, or the corrugations of boundaries are not small. 
Even in these cases, the scattering formalism is applica- 
ble, and numerical methods should prove useful. 

In the light of recent experiments on dynamical 
Casimir effect, precise computations of the effect of ge- 
ometry and motion are needed. We believe that our for- 
malism provides an efficient analytical, and possibly even 
numerical, computational tool. 
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Appendix A: Scattering matrices 

In this Appendix, we derive the scattering matrix for 
specific geometries. The Dirichlet boundary conditions 
are assumed throughout this section. 



Plate 



The perturbative scheme introduced in Sec. Ill A can 
be generalized to d dimensions. We designate the coor- 
dinates spanning the surface by xy and the normal co- 
ordinate by z. The incoming and outgoing waves are 
identified as 



in/out 



(t,x) 



1 



: exp(— zcji + ik|| • X|| =p ikj_z) , (Al) 



where ky is a (d — l)-component wavevector parallel to 

the surface of the plate, and fcj^(aj,k||) = y^w^/c^ — k|. 

Consider <I>o as the solution to the Dirichlet boundary 
problem for a static mirror. 



$o(t,x) = <i>L\,, (i,x)-Ck,,(*'^)' 



(A2) 



The scattering matrix can be computed perturbatively by 
organizing the field as $ = •I'o + $1 + • • • . The Dirichlet 
boundary condition, X|| , z + q{t, X|| )) = 0, to the first 
order is given by 

$i(t,x||,0) = -q(<,x||)a,$o(i,x||,0), (A3) 

where q{t, xy) is the boundary displacement as a function 
of time and position on the surface. Given the (time- 
dependent) value of the field <i>i on the boundary as given 
by Eq. (A3), one can compute the latter field everywhere 



in the space by using Green's theorem, 



(A4) 



where x and x' are spacetime coordinates, and Gd is a 
Green's function satisfying Dirichlet boundary conditions 
on the plate. The integral in the last equation is over a 
closed surface including x in its interior. The Green's 
function, Gd, is, in Fourier space, given by 



GD(c^,k||,z, z') 



2k ^ 



-e 



e**-^^<). (A5) 



Using Eqs. (A4) and (A5), one obtains $1, which in turn 



gives the scattering matrix as 



S'^-Hnk||+K,t^k|| = -2i q{n, K) ^/ kj_{u}, k\\) kj_{u} + n,k\\ + K) 

(A6) 
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2. Sphere 



3. Disk (cylinder in 2d) 



In spherical coordinates, the normahzed incoming and 
outgoing waves are defined as 



l1r-\[^e--^h]^-'U^)Y,„M^). (A7) 



\ c / 



The Dirichlet boundary condition for a spherical object 
in motion is 



^{t,Rf + q) = 0. 



(A8) 



In this equation, R is the radius of the sphere, f is the 
unit vector along the radius, and g is the displacement as 
a function of time and position on the sphere's surface. 
For simplicity, we assume that the sphere undergoes a 
linear (but time-dependent) motion. Equation ( |A8[ ) can 
be expanded in powers of q. To the first order, we have 

^i{t,Rf) ^ -~q{t) cos edrMt,rf)\r=R, (A9) 

where the zeroth-order solution is 

M^,^) = '^Zrn+Slii^)KTm. (AlO) 

with Si{uj) — — ^i^^^^-^/^) being the scattering matrix of a 

static sphere. The Green's function satisfying the Dirich- 
let boundary conditions on the sphere can be written as 

Gd(c.,x,x') = g E (M'H^^</c) + Si{u)h['\ujr^/c)) 



Im 



x/jp)(o.r>/c)r;™(0,0)y,:„(0',0'). 



(All) 



Green's theorem can then be applied to compute $i, or 
equivalently the scattering matrix as 



_ 2»g(») ^ 

C 



X ^/{uJ + n)uJFl (^^^ 



{uj + n)R 



(A12) 



with F defined as 



Fi{x) 



(A13) 



The constants dium are nonzero only for I' — I ± 1, 



1 Im 



'{l + m + - m + 1) 



(2/ + l)(2; + 3) 



dl^llm — dll^i, 



(AM) 



In polar coordinates, the normalized incoming and out- 
going waves are defined as (up to an irrelevant constant) 



$;"/™'(w,x) = e-'^'H^^^^\ujr/c)e'"'^. (A15) 
The boundary condition for a circular disk in motion 
is described similarly to Eq. ( A8 1 with r being a two- 
dimensional radial vector. Again we make the as- 
sumption that the object undergoes a linear (but time- 
dependent) motion along the x-axis. Equation ( |A8[ ) can 
be expanded as 



$i(i,i?f) = -q{t) cos ^drMt,rr)\r^R, (A16) 



with (j) being the angle from the x axis. The zeroth order 
solution, $0i is given by 



Mo^, x) = + SM<fZt, (A17) 



where S'm(a;) 



is the scattering matrix of 



a static disk. The Green's function subject to Dirichlet 
boundary conditions on the disk is 



Gd(u;,x,x') - ^E (i7i2)(^^</c) + 5™(c.)i7(,i)(c.r</c) 



X H^\ujry/c) e™^'^-'^') . (A18) 



With the knowledge of the field on the boundary. 



Eq. (A16), one can apply Green's theorem to obtain the 



field elsewhere in space. One then finds the scattering 
matrix 



5, 



u+nm±l,um — — — M„i I I lVlm±l I 



where M is defined as 



Mm{x) 



1 



(A19) 



(A20) 
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